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ABSTRACT 

We investigate the relation between stellar mass (MJ, star formation rate (SFR), and metallicity (Z) of galax- 
ies, so called the fundamental metallicity relation, in the galaxy sample of the Sloan Digital Sky Survey Data 
Release 7. We separate the galaxies into narrow redshift bins and compare the relation at different redshifts, and 
find statistically significant (> 99%) evolution. We test various observational effects that might cause seeming 
Z evolution, and find it difficult to explain the evolution of the relation only by the observational effects. In 
the current sample of low redshift galaxies, galaxies with different M± and SFR are sampled from different 
redshifts, and there is degeneracy between M*/SFR and redshift. Hence it is not straightforward to distinguish 
a relation between Z and SFR from a relation between Z and redshift. The separation of the intrinsic relation 
from the redshift evolution effect is a crucial issue to understand evolution of galaxies. 

Subject headings: galaxies: abundances — galaxies: evolution 



1. INTRODUCTION 

Stellar mass (M t ), star formation rate (SFR), and metal- 
licity (Z) of galaxies are essenti al parameters to under- 
stand evolution of gala x ies (e.g., iKauffmann et alj [2 003a; 
Brinchman n et all 120041: ITremonti et all l2004t ISalim et al.1 
2007). The relation b etween and Z has been studied at var- 
ious redshifts (e.g., IT remonti et al. 2004; Savaglio et al.l l2005b 
lErb et all 120061 iLiu et all 2008). Ellison et al.1 (120081) added 
the 3rd parameter to the M+-Z relation, and found the 3D cor- 
relation between M*, SFR, and Z of field star forming galax- 
ies at low redshifts (z ~ 0.1) in the Sloan Digital Sky Sur- 
vey (SDSS) galaxy sample. Recen tly, Mann ucci et alj ( 120101 
hereafter M 1 0) and lLara-Lopez et all (12010b showed that high 
redshift galaxies (z < 3) agree with the extrapolation of the 
M*-SFR-Z relation defined at low redshifts (so called the fun- 
damental metallicity relation). 

The no-evolution of the M^-SFR-Z relation with redshift 
suggests existence of a physical process which affects evolu- 
tion of galaxies in the wide range of redshift behind the re- 
lation. Some models are alread y proposed to explain the ori- 
gin of the M„-SFR-Z r elation fc>av6 et al.ll2012t iDaval etail 
l20lllYateset^l20T2T) 

Although the correlation betwe en M^, SFR, and Z is 
tested at various redshifts (M10; lLara-Lopez et al.| 
Richard et al] |20Tt1 lYabe et alJ 1201 11 iNakaiima etail 
Cresci et all 120121: IWuvts et aLl 120121) . galaxies observed at 
different redshifts typically have different M* and SFR. Hence 
the comparison of galaxy metallicities at different redshifts in 
a same range of and SFR is hardly done. Thus it is cur- 
rently difficult to distinguish a fundamental relation between 
M*, SFR, and Z that stands in wide range of redshift, from a 
series of the M+-Z relations at various redshifts. 

In this study, we separate the low redshift galaxy sample 
(z < 0.3) in which the M^-SFR-Z relation is defined into nar- 
row redshift bins, and investigate redshift evolution of the re- 
lation. Throughout this paper, we assume the fiducial cosmol- 
ogy with n A = 0.7, Q m = 0.3, and H = 70 km s" 1 Mpc" 1 . 
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2. THE GALAXY SAMPLE 

We draw galaxies for our analysis from the SDSS MPA- 
JHU Data Release 7 catalog^ (hereafter the MPA-JHU cat- 
alog), with similar selection criteria to that in M10 as fol- 
lows. We select galaxies (1) at redshifts between 0.07 and 
0.3, (2) the signal-to-noise ratio (S/N) of Ha line > 25, (3) 
A v < 2.5, (4) Ha to H/3 flux ratio (Ha/H/3) > 2.5, and (5) 
not str ongly affected by active galactic nuclei (AGN) accord- 
ing to Kauffmann et al. (2003b). These selection criteria are 
identical to those in M10. 

The target galaxy selection method for the spectroscopic 
observation in the SDSS provides highly uniform and com- 
plete sample of galaxies above applied magnitud e limit that 
is r-b and Petrosian magnitude m,.p et < 17.77 ( Strau ss"et al.l 
2002). However, galaxies fainter than this limit which have 
been selected by dif ferent selection met hods (e.g. luminous 
red galaxy sample. lEisenstein et al.ll200l1) are also included in 
the MPA-JHU catalog. To avoid possible selection effect in 
those galaxies, we impose following selection criteria in ad- 
dition to the M10 criteria. Galaxies in our sample (6) have 
m r ,Pet < 17.77, (7) are selected for spectr oscopic targets as 
galaxy candidates (Stoughton et al. 2002), and (8) spectro- 
scopically confirmed as galaxies by the SDSS spectroscopic 
pipeline (see the SDSS DR7 websitsfj for detail). Finally 
~ 110000 galaxies are left in our sample. We note that ~ 
90% of the galaxies in our sample is at z < 0.15, although we 
collect galaxies at redshifts up to 0.3. 

We us e M* and SFR of galaxi e s listed in the MPA-JHU 
catalog ([K auffman neTaT] l2003at iBrinchmann et alJ l2004t 
ISalim et al.l2007ir and me asure metallicity of the galaxies us- 
ing empirical method by Maiolin o et al.l (120081) . Following 
M10, we use R 2i = ([OII]A3727+[OIII]A4958,5007)/H/3, and 
[NII]A6584/Ha as the indicators of metallicity. The line flux 
ratios are corrected for extinctions us i ng Ha /H/3 and the ex- 
tinction curve given by ICardelh et al.l ([1989). When both in- 
dicators are applicable to a galaxy, we use mean of the 12 + 
log(0/H) values given by the two indicators, as far as the two 
values agree with each other within 0.25 dex. When the two 

2 http://www.mpa-garching.mpg.de/SDSS/DR7/ 

3 http://www.sdss.org/dr7/ 
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indicators disagree, we remove the galaxy from our sample. 
Only 3% of galaxies in our sample is rejected here, similarly 
to the case in M10. 

3. CORRELATION OF SFR AND METALLICITY 

We show the correlation between SFR and metallicity of the 
galaxies in mass ranges log M^/Mq = 10.0 ±0.05, 10.25 ± 
0.05, 10.5 ±0.05, 10.75 ±0.05, and 11.0 ±0.05 in figure Q] 
We separate the galaxies into SFR bins with Alog SFR = 0.1 
dex, and measure median of metallicities in each bin. The cor- 
relation between SFR and metallicity is not clear among the 
galaxies with log M*/M > 10.5 (the top three panels). On 
the other hand, among the galaxies with log M+/Mq < 10.25, 
higher SFR galaxies have lower metallicity than lower SFR 
galaxies (the bottom two panels), consistently to the results of 
M10. 

We also separate the galaxies into redshift bins with Az = 
0.02, and compare the SFR-Z relations at different redshifts 
(figure [TJ. The plotted redshift bins are different in different 
mass ranges, due to the difference of redshift distribution of 
the galaxies in each mass range (the top panel of figure 
Each data point in figure Q] contains more than 50 galaxies. 
Although systematic difference of the SFR-Z relation at dif- 
ferent redshifts is not clear in the mass ranges log M+/Mq > 
10.75, the galaxies with with log M^/Mq < 10.5 show red- 
shift evolution of the SFR-Z relation. 

In the mass ranges log M+/M Q < 10.5, galaxies with higher 
SFR shows larger chemical evolution. Making the negative 
correlation between SFR and Z steeper at higher redshift. The 
low redshift sample of galaxies with log M^/Mq = 10.5 shows 
a positive correl ation betwee n SFR and Z which is consistent 
to the results of lYates et aD (120 121) . while the high redshift 
galaxies with the same mass range shows the negative corre- 
lation similarly to the case in lower mass ra nges. 

As discussed in Brisbin ^fe Harwij (120121) . the extent of the 
metallicity evolution with redshift shown in figure [T]is small 
(< 0.05 dex) due to the narrow redshift range we can investi- 
gate with our sample. However, most datapoints which show 
the redshift evolution in figure[T]have the error of mean < 0.01 
dex indicating that the evolution is statistically significant. We 
discuss stat istic al significance of the evolution of the relation 
further in § |4~T1 

It should be noted that a combination effect of fiber cover- 
ing fraction and metallicity gradient in a galaxy may cause 
seeming metallicity evolution with redshift. In the SDSS, 
galaxy spectra are obtained only in the 3 arcsec fiber aperture. 
The fiber aperture can contain larger area of a target galaxy 
at higher redshifts, while galaxies tend to hav e lower metal- 
licity in their outskirt than at their cent er (e.g., Zaritsky et al. 
1994). We discuss this issue in § 14.21 There are also other 
observational effects that may affect the redshift evolution of 
the M*-SFR-Z relation. The limiting magnitude of the SDSS 
spectroscopic target selection m r ,p et < 17.77 and/or the Ha 
S/N > 25 threshold in our sample selection may cause some 
redshift dependent sampling bias. The larger noise in galaxy 
spectra at higher redshift is also a possible sou rce o f artif icial 
redshift effect. We discuss these effects in § 14.31 and § 14.41 
respectively. 

In figure Q] the galaxies with log M^/Mq > 10.5 and 
log SFR [M Q yr _1 ] < 0.5 show possible negative evolution of 
metallicity (higher metallicity at higher redshift), although the 
error of mean is large (the sample size is small). One pos- 
sible explanation for the negative chemical evolution of the 
low SFR galaxies is decrease of SFR in galaxies with high 
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FIG. 1. — The correlation between SFR and Z of the galaxies with 
log M*/M Q = 11.0 ±0.05, 10.75 ±0.05, 10.5 ±0.05, 10.25 ±0.05, and 
10.0 ±0.05 from top to bottom, respectively. The data points connected with 
solid lines (black) represents the median metallicity of the galaxies in the 
whole range of redshift 0.07 < z < 0.3 in each SFR bin. The data points con- 
nected with dashed (blue) and dotted (red) lines represent the median metal- 
licities of the low and high redshift samples. The vertical error bars represent 
the error of mean of the metallicity distribution in each bin (the scatter of 
the sample divided by the square root of the sample number). The horizon- 
tal error bars represent the width of the SFR bin. The low and high redshift 
data points are slightly shifted sideway for visibility. The grey dots represent 
individual galaxies. A colored version of the figure is available in the online 
journal. 
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FIG. 2. — Top panel: the redshift distribution of the galaxies in our sample 
with stellar mass: log M*/Af© = 10.0 ± 0.05, 10.5 ± 0.05, and 11.0 ± 0.05. 
Bottom panel: the redshift distributions of the log M*/Mq = 10.5 galaxies, 
with log SFR [M yr' ] = 0.2 ± 0.05 and 0.8 ± 0.05. 

SFR and low metallicity. When SFR and metallicity are neg- 
atively correlated, some of high-SFR low-metallicity galax- 
ies may decrease their SFR before significant chemical evo- 
lution decreasing the median metallicity of low SFR galaxies, 
while others continue star formation undergoing chemical en- 
richment and increasing the median metallicity of high SFR 
galaxies. 

To investigate the redshift evolution further, we plot the me- 
dian metallicities of the galaxies with log M+/Mq = 10.5 ± 
0.05 and various SFR as functions of redshift in figure [3] 
The metallicities of the high SFR galaxies (log SFR [M Q yr _1 ] 
> 0.5) is lower at higher redshifts, while the metallicities of 
the low SFR galaxies shows negative or no evolution, con- 
firming the results in figure Q] 

In the bottom panel of figure [2] we show redshift distribu- 
tions of the galaxies with log = 10.5, and log SFR 
[M Q yr _1 ] = 0.2 & 0.8. The galaxies with different SFR are 
sampled from different redshifts, as well as in the case of dif- 
ferent shown in the top panel. In figureQ] the SFR-Z rela- 
tions for the whole redshift range (0.07 < z < 0.3) follows the 
low redshift relation in the low SFR end, while metallicity of 
the high SFR galaxies follows the high redshift relation, con- 
sistently to the difference of the redshift distributions between 
the low and high SFR galaxies. 

4. DISCUSSION 
4.1. Kolmogorov-Smirnov Test 

To examine the significance of the redshift evolution we 
find in §[3] we perform the two-sample Kolmogorov-Smimov 
(KS) test, which calculates the probability that two samples 
can be drawn from a same probability distribution function. 
In figure [4] we show metallicity distributions of the galaxies 
with the four sets of and SFR, which show notable redshift 
evolution in figureQ] at two different redshifts each. 

The positive metallicity evolution (lower-metallicity at 
higher-redshift) in the mass ranges log M+/Mq < 10.5 is sta- 
tistically significant to high confidence level > 99% (the top 
3 panels of figure 0J. The statistical significance of the evo- 
lution of the galaxies with log M+/Mq ~ 10.25 and 10.5 is 
very high with the KS test probability P K s = 7.7 x 10" 13 and 
2.1 xlO" 7 . 

Although the statistical significance of the evolution of the 
galaxies with log M+/Mq = 10.0 is lower than in the cases 
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FIG. 3. — The metallicity evolution of the galaxies with log M^/Mq = 
10.5 ± 0.05. The data points represent the median, the la scatter (thin error 
bars), and the error of mean (thick error bars) of the galaxy metallicities. The 
width of each bin is 0.02 dex and 0.1 dex for redshift and SFR, respectively. 
The data points are slightly shifted sideway for visibility. 

of the galaxies with log M±/M Q ~ 10.25 and 10.5 due to 
the very narrow range of redshift where we can investigate 
log M*/M© = 10.0 galaxies, the evolution in this mass range is 
also larger than statistical error {Pks = 1 .5 x 10~ 3 ). The extent 
of the evolution of the galaxies with log M+/Mq = 10.0 is sim- 
ilar to that of the log M*/M = 10.5 galaxies (~ 0.03 dex), al- 
though the redshift range investigated is 3 times smaller. This 
suggests that less massive galaxies are more rapidly evolving 
in metallici ty, consistently to the findings of previous stud- 
ies (e .g., ISavaglio et al.l 120051 but see also Lamareill e et alj 
12001 . 

The possible negative metallicity evolution of the high- 
mass, low-SFR galaxies is statistically less significant than 
the positive evolution of the high SFR galaxies (the bottom 
panel of figureHJi, as expected from the large error of mean in 
figureQ] 

4.2. The Effect of Fiber Covering Fraction 

The combination effect of fiber covering fraction and metal- 
licity gradient in a galaxy may cause seeming metallicity evo- 
lution with redshift. The SDSS spectroscopic fiber (3 arcsec) 
contains larger area of a target galaxy at higher redshift, and 
galaxies tend to have lower metallici ty in their outskirt than at 
their center (e.g., Zarits kvet al.lll994B . 

We show the relation between the fiber covering fraction 
(the fiber to total flux ratio in r-band) and the derived metal- 
licity of the galaxies with log M*/M = 10.5 and log SFR 
[M Q yr -1 ] = 0.8 in figure [5] The galaxies with larger covering 
fraction have higher-metallicity, contrary to what is expected 
from the metallicity gradient. This trend may be a result of 
the negativ e correlation betwee n galaxy radius and metallic- 
ity found bv lEllison et al.l (120081) . The absence of the expected 
negative correlation between the fiber covering fraction and 
the metallicity suggests that the fiber covering fraction is de- 
termined by intrinsic size of galaxies rather than their red shift. 
We note that the galaxies have significant variation of radus 
(6-15 kpc, Petrosian) even in the narrow bin of and SFR. 

It should be noted that the gradient effect may be hidden in 
figure [5] canceled by the effect of the radius-metallicity cor- 
relation, and hence the absence of the gradient effect in the 
plot doesn't completely rule out the effect. However, if the 
redshift evolution of metallicity results from the effect of the 
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FIG. 4. — The cumulative metallicity distributions of the galaxies with the 
four sets of M* and SFR, at low and high redshifts. The top 3 panels show the 
high SFR galaxies whose metallicity is lower at higher redshift, while the bot- 
tom panel shows the low SFR galaxies whose metallicity is (possibly) higher 
at higher redshift. In each panel, solid and dashed histograms represent low 
and high redshift samples, respectively. The plotted redshifts are different for 
different set of M* and SFR, due to the difference of the redshift distributions 
of the galaxies with different M* and SFR (see figure[2). 

metallicity gradient, we expect larger evolution for galaxies 
with larger radius and the gradient. Considering that galaxies 
with larger M t tend to have larger radius, and dwar f galaxies 
have flatter metallicity gradient than spirals (e.g.. iLee et alj 
2006), the finding of the metallicity evolution only in high 
specific SFR galaxies is opposite to the expected trend in the 
case of the gradient effect. Hence we consider it is difficult to 
explain the metallicity evolution found in this study only by 
the effect of the metallicity gradient. 

4.3. Continuum and Ha Line Luminosity Limits 

As mentioned in §|2] the SDSS spectroscopic target galax- 
ies are magnitude limited m r p et < 17.77, and we select galax- 
ies with Ha S/N > 25 for this study. Hence the broad band 
limiting magnitude and/or the Ha flux limit may cause some 
redshift dependent sampling bias. 

In figure [6] we plot M r and Ha luminosity distributions of 
the galaxies with log M*/M Q = 10.5 and log SFR [M Q yr _1 ] 
= 0.8 at redshift 0.08 and 0.14. With the fixed and SFR, 
the galaxies at z = 0.14 are clearly brighter than the z = 0.08 
galaxies in r-band, suggesting that the r-band limiting magni- 
tude plays an important role in this sample. On the other hand, 
the Ha luminosity distributions at the two different redshifts 
are not largely different. We note that galaxies with log SFR 
[Meyr 1 ] = 0.8 typically have Ha S/N - 60 at z = 0.14, and 



o 

a 9.02 



-0.80 -0.60 -0.40 

lo 8 F nber /F tota | 

FIG. 5. — The median, the 1<t scatter, and the error of mean of the galaxy 
metallicities with log M^/Mq = 10.5 and log SFR [MQyr -1 ] = 0.8, as a func- 
tion of the fiber covering fraction. Here the fiber covering fraction is the ratio 
of flux in fiber to total flux of the object measured in r-band. The width of 
the fiber covering fraction bins is 0. 1 dex. 

hence few galaxies in this SFR range are rejected by the Ha 
S/N limit by the Ha S/N threshold. 

To discuss the effect of the limiting magnitude, we plot dis- 
tributions of the galaxies with log M^/Mq = 10.5 in some red- 
shift bins on the SFR versus r-band mass-to-luminosity ratio 
(M*/L„ r ) plane and the metallicity versus M t /L„. r plane in 
figure 

In the left panel of figure[7] M+/L v>r of the galaxies clearly 
correlates with SFR, causing the difference of redshift distri- 
bution between the galaxy samples with different SFR (the 
bottom panel of figure |2). On the other hand, no clear cor- 
relation between metallicity and M*/L U j is seen in the right 
panel. Furthermore, it is also notable that galaxies at higher 
redshifts have systematically lower metallicity than lower red- 
shift galaxies with similar M^jL VJ . 

To test the effect of the limiting magnitude further, we 
separate the galaxies with log M+/Mq = 10.5 and log SFR 
[M Q yr _1 ] = 0.8 at redshift 0.08 into two subsamples. One is 
bright (M r < -21.5), and the other is faint (M r > 21.5). We 
compare the metallicity distributions of the bright/faint sam- 
ples at z= 0.08 and the sample at z = 0.14 in figure [8] The 
galaxies in the bright sample at z = 0.08 have similar M r to 
the galaxies at z = 0.14. Although the bright sample is small, 
the KS test between the bright sample and the z = 0. 14 sample 
indicates the metallicity distributions of the two sample is sig- 
nificantly different (Pks = 0.037), while the faint sample and 
the bright sample is broadly consistent (Pks = 0.31). Hence 
it is unlikely that the limiting magnitude effect is the primary 
source of the redshift evolution of the M^-SFR-Z relation. 

4.4. Noise Effect and Mean Spectra 

With similar and SFR, spectra of high redshift galax- 
ies would have lower S/N than those of low redshift galaxies. 
This expected trend possibly affect the metallicity estimate in 
a redshift dependent way. We also note that figure [6]suggests 
that high redshift galaxies have smaller Ha equivalent width 
than low redshift galaxies, and hence it is possible that high 
redshift galaxies suffer more from uncertainties of the con- 
tinuum subtraction than low redshift galaxies. To reduce the 
noise effect, we compose mean spectra of the galaxies with 
log MJM Q = 10.5 ±0.05 and log SFR [A^yr" 1 ] = 0.8 ±0.05 
at redshifts 0.08 ±0.01, 0.10±0.01, 0.12±0.01, 0.14±0.01, 
and 0.16 ± 0.01, and measure [NII]A6584/Ha and R 2 i of the 
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mean spectrum at each redshift. 

We average spectra collected from the SDSS data archive 
serveiQ after shifting wavelength to the center of each red- 
shift bin. To reduce the effect of stellar absorption lines on 
the emission line measurement, we fit continuum of the mean 
spectra with stellar spectral energy distribution (SED) mod- 
els, and subtract the stellar component models form the mean 
spectra. We perform the SED fitting with the SEDfit software 
package (ISawickill2012t [Sawicki & Yedll998l) which utiliz e 
the population synthesis models o flBruzual & Cha riot (2003). 
We examine five cases of star formation history: simple stellar 
population, constant star formation, and exponentially decay- 
ing star formation with r = 0.2, 1.0, and 5.0 Gyr assuming 
stellar metallicity to be Z Q . 

Two mean spectra and their best fit stellar SED models are 
shown in the left panel of figure [9] The mean spectra achieves 
very high S/N (typically ~ 100 in the highest redshift bin). 
The best fit parameters of the stellar SED models indicate 
smaller M„ and SFR than log M*/M = 10.5 and log SFR 
[M yr _1 ] = 0.8, but we note that the spectra only represents 
stellar populations in the spectroscopic fiber. In the right panel 
of figure [9] we show a close up view of high order Balmer 
lines in the mean spectra. Although the center of each ab- 
sorption line is affected by the corresponding emission line, 
the stellar SED models reproduce Balmer absorption lines in 
their wings. 

After subtracting the stellar SED models, we fit each line 
with linear continuum plus single gaussian to measure flux of 
each line (we fit [OII]A3727 doublet with linear continuum 
plus double gaussian). We show the line fittings at z = 0.08 
and 0.14 in figure [10] and [TT] When normalized to the peaks 
of the gaussian models of Ha, the spectrum at z = 0.08 has 
slightly stronger [NII]A6584 than the spectrum at z = 0.14 
(figure [TUll. We also fit [NII]A6548 to avoid that the line 

4 http://das.sdss.org/www/html/ 



affects the estimate of the residual continuum, although we 
don't use the line as a metallicity indicator. One can also 
find that [OII]A3727 is slightly weaker at z = 0.08 than at 
z = 0.14 (figure [UJ, contrary to the case of [NII]A6584. We 
note that the continuum level underneath [OII]A3727 may be 
affected by high order Balmer emission lines whose peaks are 
not resolved in the mean spectra, and it may affect R23, while 
[NII]/Ha is free from this effect. 

The line fluxes obtained by the gaussian fit are listed in ta- 
ble [TJafter_corrected for extinction using Ha/H/3 ratio and the 
Carde lli et al.l (|T989) extinction curve. The redshift evolution 
of the metal indicating line ratios is plotted in figure [12] to- 
gether with the median line ratios of the galaxies with the 
same M*, SFR, and redshift to those of the mean spectra. Note 
that the line fluxes of the mean spectra are obtained by our 
gaussian fit, while the median line ratios are obtained using 
the line fluxes listed in the MPA-JHU catalog. 

Galaxies at higher redshifts show smaller [NII]/Ha and 
larger R23, in their mean spectra. The relations between 
the line ratios and Z are ^log([NII1/Ha)/<jl ogZ = 0.81 and 
Qflogj ? M MogZ = -1.5 at 12 + log(0/H) = 9.0 dMaiolino et alj 
2008), and hence the evolution of the line ratios in figure [T2l 
is broadly consistent to the metallicity evolution found in §[3] 
The median line ratios also evolve with similar manner to the 
line ratios of the mean spectra, although they indicate smaller 
metal-to-hydrogen line ratio than the results with the mean 
spectrum, possibly due to the difference of spectrum fitting 
method. Thus we conclude that the metallicity evolution is 
not a result of noise effect. 

5. SUMMARY 

We have investigated the evolution of the M*-SFR-Z re- 
lation at z < 0.3. Although the redshift range we investi- 
gated is narrow, we found metallicity of the high SFR galaxies 
(log SFR [Moyr" 1 ] > 0.5) with log M„/M < 10.5 is evolv- 
ing with > 99% statistical significance. We have examined 
the observational effects which may cause seeming evolution 
of metallicity: the fiber aperture effect, the sampling bias by 
the limiting magnitude, and the noise in the spectra. We found 
it is difficult to explain the evolution of the M^-SFR-Z rela- 
tion only by the observational effects, although some effects 
are not completely ruled out. 

In the current galaxy sample at low redshifts, galaxies with 
different and SFR are sampled with different redshift dis- 
tributions, and there is a degeneracy between MJSFR, and 
redshift. Hence it is difficult to clearly separate SFR depen- 
dence of metallicity from redshift dependence. The metallic- 
ity evolution of galaxies with same and SFR at z < 0.3 
found in this study suggests that the currently known M + - 
SFR-Z relation arises from a combination of the intrinsic re- 
lation and the redshift evolutio n effect. Previous studies (e.g., 
M10: ILara-Lopez et al.l 12010) showed that z > 0.5 galaxies 
agree with the extrapolation of the relation defined by z < 0.3 
galaxies, suggesting that the relation is not evolving with red- 
shift. However, if the redshift evolution effect is included in 
the relation at z < 0.3, the agreements of high redshift galaxies 
don't necessarily mean that the relation is not evolving. The 
agreements may result from the extrapolation of the redshift 
evolution effect. 

We need to investigate metallicity distribution of galaxies 
with same and SFR at different redshifts, to distinguish the 
intrinsic relation from the redshift evolution effect. However, 
it is currently possible only for narrow range of M + , SFR, and 
redshift. The extent of the evolution we find in the current 
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FIG. 7. — Left panel: the distribution of the galaxies with log M*/Mq = 10.5, on the SFR versus mass-to-luminosity ratio plane. The distributions of 
the galaxies at z = 0.08,0.12 and 0.16 ±0.01 are shown with solid (blue), dashed (magenta), and dotted (red) contours, respectively. The contours indicates 
rf 2 ngai/rflogSFRrflog(M*/L„, r ) = 1,3,9... [dex~ 2 ]. The r-band solar luminosity is L vj ,q = 5.9 x 10 18 erg sec" 1 Hz" 1 dBlanton et alj|2003l) . Right panel: same 
as the left panel, but on the metallicity versus mass-to-luminosity ratio plane. 



TABLE 1 

The Emission Lines in the Mean Spectra 



redshift [OII]A3727 H/3 [OIII]A4959 [OIII]A5007 Ha [NII]A6584 



0.08 


2.72 X 10~ 


14 


1.50 X 10" 


14 


1.50 X 10- 


15 


4.15 X 10- 


15 


3.86 x 10" 


14 


1.51 x 10" 


14 


0.10 


1.91 x 10- 


14 


9.98 x 10" 


15 


9.95 x 10- 


16 


2.88 X 10- 


15 


2.57 x 10" 


14 


9.80 x 10" 


15 


0.12 


1.34 x 10' 


14 


6.86 x 10- 


15 


7.32 x 10- 


16 


2.10 X 10- 


-15 


1.78 x 10" 


14 


6.40 x 10" 


15 


0.14 


9.18 x 10" 


15 


4.76 x 10- 


15 


5.38 x 10- 


16 


1.46 x 10" 


15 


1.24 x 10" 


14 


4.46 x 10- 


15 


0.16 


8.14 x 10" 


15 


3.58 x 10- 


15 


4.37 x 10- 


16 


1.27 X 10- 


-15 


9.36 x 10" 


15 


3.31 x 10- 


15 



NOTE. — The extinction corrected line fluxes of the mean spectra of the galaxies with log M^/Mq 
10.5 ±0.05 and log SFR [M yr'] = 0.8 ± 0.05. The units are [erg s" 1 cm" 2 ]. 




FIG. 8. — The cumulative metallicity distributions of the galaxies with 
log M^/Mq = 10.5 and log SFR [M yr'] = 0.8, at the two different red- 
shifts. The low redshift sample is separated into the bright (M r < —21.5) and 
the faint (M r > 21.5) subsamples. The galaxies in the bright sample have 
similar M r to that of the high redshift sample, where faint galaxies are not 
selected as spectroscopic targets. The KS test probabilities between the dis- 
tributions are indicated together. 



sample is very small due to the very narrow range of redshift 



we can investigate. Hence it is difficult to robustly exclude 
possibility of systematic effects in our results. 

Deeper wide field spectroscopic survey which covers wider 
range of redshift is necessary to test the evolution effect more 
robustly. For example, to study similar range of and SFR 
at z ~ 1 with a comparable sample number (or a survey vol- 
ume, > 10%) to this study, the limiting magnitude must be 
5 magnitude deeper than the SDSS sample (m < 23), and the 
survey field must be > 20 deg 2 (0.3% of the SDSS DR7 field). 
Note that the comoving volume element is ~ 37 times larger 
at z = 1 than at z = 0.1. Future spectroscopic surveys, such 
as Subaru Measurement of Images and Redshifts (SuMIRe) 
project, will achieve such requirements. At redshifts z < 1, the 
limiting magnitude can be brighter but the survey field must 
be wider. The robust separation of the intrinsic M^-SFR-Z 
relation from its evolution will open a new way to understand 
galaxy evolution. 



We are grateful to the MPA/JHU group for making their 
galaxy catalog public. We thank K. Yabe for helpful dis- 
cussions. Thanks are also due to our referee, M. A. Strauss, 
whose thoughtful comments largely improved this paper. 
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FIG. 1 1 . — Same as figure fTol but for R23 lines. The right vertical axis indicates f\ normalized to the peak of H/3 line. 
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FIG. 12. — The redshift evolution of the metal indicating line ratios of the 
mean spectra. The median line ratios of the galaxies with the same M», SFR, 
and redshift are plotted together. Each line is detected with S/N = 100-1000 
in the mean spectra, where the vertical error bar in the left corresponds to S/N 
= 100. The vertical error bars of the median line ratios indicate the error of 



